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A microscopic theory for the spin triplet Cooper pairing in non-centrosymmetric supercon- 
ductors like CePtaSi and CeTSi3 (T=Rh, Ir) is presented. The lack of inversion symmetry 
leads to new anomalous spin fluctuations which stabilize the triplet part in addition to the 
singlet part originating from the centrosymmetric spin fluctuations. It is shown that both parts 
have similar nontrivial momentum dependence of Ai type. Therefore the mixed singlet-triplet 
gap function has accidental line nodes on both Fermi surface sheets which are stable as func- 
tion of temperature. This gap function explains the salient features of CePtaSi and CeTSi3 
superconductors. 
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Since its discovery in the heavy fermion compound 
CeCu2Si2, 1 unconventional superconductivity has at- 
tracted much attention in the field of correlated elec- 
tron systems. For systems with inversion symmetry, the 
pairing states of superconductivity are classified as spin- 
singlet of even-parity or spin-triplet of odd-parity. Re- 
calling that a phonon-mediated mechanism is established 
for s-wave spin-singlet superconductivity, the possibility 
of spin-triplet superconductivity, as realized, e.g. in UPt3 
and Sr2Ru04, is considered as a characteristic feature of 
the correlated electron system. 

It has been suggested that the spin-triplet supercon- 
ductivity is destroyed by breaking the inversion symme- 
try. 2 However, recently, the unconventional superconduc- 
tivity has been observed in several non-centrosymmetric 
heavy fermion compounds like CePtaSi, 3 CeRhSi3, 4 
CeIrSi3, 5 and CeCoGe3. 6 These compounds have tetrag- 
onal crystal structures without the inversion symmetry, 
and the superconductivity is found around each antifcr- 
romagnetic phase. 4-7 The evidence of unconventional su- 
perconductivity is at least shown by the existence of a 
line node of the superconducting gap in CePtaSi, which is 
observed by NMR, 8 thermal transport, 9 and penetration 
depth. 10 Furthermore, the common remarkable property 
of the unconventional superconductivity is the large up- 
per critical magnetic field exceeding the Pauli limiting 
field. 3 ~ 5 ' n The observed upper critical field has opened 
up the possibility of the spin-triplet superconductivity in 
non-centrosymmetric systems in contrast to the theoret- 
ical expectation. 2 

The lack of inversion symmetry of non- 
centrosymmetric systems leads to an antisymmetric 
spin-orbit interaction. 12 A general superconducting or- 
der parameter is described by the even spin-singlet gap 
function ip(k) and odd spin-triplet gap function d(k). 13 
The coexistence of ip(k) with d(k) is unavoidable in 
principle in the non-centrosymmetric superconductor 
because of the lack of inversion symmetry. For the 
non-centrosymmetric superconductor with a large 



upper critical field, Sigrist et al. have reconsidered the 
possibility of spin-triplet pairing and concluded that a 
Cooper pair satisfying d(k) || gk (gk: Rashba field) does 
not suffer the suppression of pairing coming from the 
lack of inversion symmetry. 14 Based on this spin-triplet 
state, many theoretical proposals have been made for 
the superconducting properties. However, the pairing 
mechanism remains an unsettled problem except for a 
few attempts. 15-17 

On the other hand, a spin fluctuation property char- 
acteristic of the non-centrosymmetric system has been 
developed recently. 18 Unlike the centrosymmetric sys- 
tem, the anomalous spin fluctuations do not vanish in 
the non-centrosymmetric system under the condition, 
that the symmetry of the momentum dependence is 
equivalent to the symmetry of spin-product included 
in the spin fluctuation. Furthermore, not only usual 
(centrosymmetric) spin fluctuations but also anomalous 
(non-centrosymmetric) spin fluctuations enhance on ap- 
proaching the magnetic instability. Therefore, it is ex- 
pected that the latter play an essential role in the pairing 
mechanism of superconductivity. 

In this Letter, we study the superconducting transi- 
tion in a correlated electron system without inversion 
symmetry. The roles of usual and anomalous spin fluctu- 
ations in the superconductivity and the nodal structure 
of gap functions are examined within the weak-coupling 
theory. We also compare the superconductivity of cen- 
trosymmetric systems with that of non-centrosymmetric 
ones around the magnetic instability. 

The latter is described by the following Hamiltonian 
H = Hq + Hi for correlated electrons: 

',(0.1) 



Hp = ^ [(£k - n)vo + gk • o]<T<T>c 



Ckt 



Hi 



kcrcr' 



(0.2) 



where Ck CT and ct are annihilation and creation opera- 
tors of an electron with a momentum k and a spin a. 
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Here, £k and \i are the energy dispersion of electrons 
and the chemical potential, respectively, while gk= g-k 
describes the Rashba field coming from the antisym- 
metric spin-orbit interaction, which breaks the inver- 
sion symmetry. Then, eigenenergies of Hq are given by 
£k± = £k ± |gk| — [i- In Hi, U is the screened on-site 
interaction. 

In the following, we consider a two-dimensional tetrag- 
onal system with a dispersion energy £k — 2t\(cosk x + 
cos k y ) + 4<2 cos k x cos k y and a Rashba-field gk = 
3 (sin k v , — sink x , 0), which is a periodic form of the 
simplest one (k y , -k x , 0). 14 Chosing parameters as 
i2/ti=0.35 and g/t\=0.2, we can reproduce quasi-two 
dimensional Fermi surfaces of CePtaSi obtained by band 
calculations. 19,20 Furthermore, the reason why the two 
dimensional system is used instead of the realistic three 
dimensional system is based on the following knowl- 
edge. It is known that the unconventional superconduct- 
ing phase shrinks with increasing the dimensionality in 
the centrosymmetric system. 21 Recalling that the su- 
perconducting transition temperature 1.6 K in a non- 
centrosymmetric system CeIrSi3 is in the same order as 
the highest superconducting transition temperature 2.7 
K of CeCoIn 5 among Ce-based compounds, 22 it is sug- 
gestive that unconvensional superconductivity is more 
favorable in two dimensional system than in three di- 
mensional one, independent of the inversion symmetry. 

From the specific heat data 23 and 2A/T C value used 
to fitting to the NMR relaxation rate 24 in CeIrSi3, the 
strong-coupling theory is suggested. However, the weak- 
coupling theory will be sufficient to examine the mech- 
anism of superconductivity. In general, the matrix gap 
function A(k) is decomposed into a spin-singlet tl>(k) and 
a spin-triplet d(k)-vector as A(k) = [i/j(k)ao + d(k) • 
a]ia y . 13. In non-centrosymmctric superconductors, only 
spin-triplet component satisfying |d(k) -gk| =|d(k)||gk| 
is not affected by the suppression of pairing coming from 
the lack of inversion symmetry. 14 Then, the spin-triplet 
component will be given by d(k)=^(k)gk, where the 
symmetry of momentum dependence of 4>(k) is the same 
as that of the spin-singlet V(k). 16 ' 25 Under this condi- 
tion for d(k), the normal and anomalous matrix Green's 
functions, G(k, \u n ) and F(k, iu> n ), respectively, are de- 
fined with gk=gk/|gk| as 

G(k,iw„) = G + (k,kj„)o- + G_(k,kj„)g k • a ,(0.3) 

F(k, iw„) = [F + (k, iw„)<7 + F_(k, iw„)g k • 

G±(k, iw n ) and F±(k, iu) n ) are given by 



G±(k,iw„) = - 



F±(k,iu>„) = - 



£k+ 



-lCJr, 



£k- 



El 



k+ 



^k+ 



± 



A k 



(0 T 5) 
(0.6) 



A 2 



with Ak±='0(k) ± </>(k)|g k and E^-y o k± 

Within the weak-coupling theory for superconductiv- 
ity, only the static spin susceptibility is required. For 
simplicity, we neglect the feedback effect of supercon- 
ductivity on the spin fluctuation below T c . Therefore, 
in this case the spin fluctuation is just the static spin 



susceptibility in the normal state. Including the electron 
repulsion U within RPA, the matrix of static spin sus- 
ceptibility is calculated with the above Green's function 
as x(q) = [l — 2Ux^°\q)] X^H*!)) wnere the matrix 
element X^(q) for a, /3=x, y, z is given by 

_L \p r af3 (\c- f( £k ^ ~ ^( £k +q«) 

8iV 



xS(q) 



-, (0.7) 



£k+q{ - £k( 

with the Fermi distribution function f(e) and a vertex 

(k;q)=<5 a>/3 (l-£Cgk-gk +q ) (0.8) 

+£,( (SkaSk+q/J + <?k/3ffk+qa) — £a/3 7 i (£ffk+q 7 - Cffk 7 ) ■ 

Similarly, the charge fluctuation is described by Xcc(q) = 
x£ ) (q)/[l+2J7xc?(q)] withxi?(q), which is obtained by 
the replacement of T^(k; q) by T^(k; q)= l+£Cgk'gk+ q 

(0)/ \ 

in Xa/jW- 

We comment on the spin fluctuations in the tetragonal 
system. In the centrosymmetric system, the matrix x(q) 
is diagonal. We denote Xzz{q) and (Xxx(q) + Xw(q))/ 2 
as usual (centrosymmetric) spin fluctuations. On the 
other hand, in the non-centrosymmetric case, (Xxx(q) — 
Xw(q))/2, (Xx V (q) + Xj/x(q))/2, (x V z(q) - Xz V (q))/2i, 
and (xzx (q) — Xxz (q) ) /2i also remain non-zero with char- 
acteristic q-dependences of <? 2 — q 2 -, q x q y -, q y -, and q x - 
type, respectively. 18 We call these four contributions 
as anomalous spin fluctuations. The characteristic q- 
dependences of anomalous (non-centrosymmetric) spin 
fluctuations are caused by the Rashba field gk under the 
symmetry constraint. Considering the matrix x(l)i n °t 
only usual spin fluctuations but also anomalous spin fluc- 
tuations develop around a magnetic instability. 

We now examine the effect of both usual and anoma- 
lous spin fluctuations on the pairing mechanism in the 
non-centrosymmetric tetragonal system. Using the stan- 
dard procedure, 26-29 the following gap equation is ob- 
tained, 



[ ^(k) d X (k) dy(k) ] 



N ^ 



V xs (q) 

Vysid) 



V sx {q) 
V xx {q) 



Vx V {q) 

v yy{ 



Vv V (q) 



F s (k-q) 
F x (k-q\0.p) 
F y (k-q) 



where V^(q) is the pairing interaction due to corre- 
sponding fluctuation exchange given by 

-U 2 [xzz(q) + Xxx(q) + Xw(q) - Xcc(q)] - u 

^ 2 [Xcc(q) + Xzz(q) - {Xxx(q) - X TO (q)}], 

c /2 [xcc(q) + Xzz(q) + {xxx(q) - x TO (q)}], 

-U 2 [Xxy(ci) + Xyx(ci)}, (0.10) 



V ss (q) 

Vxx (q) 

Vx V {q) 
V sx (q) 
V sy (ci) 

Similar relations between the pairing interactions and 
spin fluctuations in non-centrosymmetric systems are ob- 
taind from a different approach. 30 In order to get a favor- 
able form of the gap equation, F s (k) and F a (k) are intro- 
duced by T J2 n E(k, iui n ) = [F s (k)d- +F(k)-a}iay, where 
F s (k)=[V(k)<Mk) + d(k) • g k ^-(k)]/2 and F a (k)= 



V yx (q) 



-V xs (q) = iU 2 [x yz {q) ~ Xzy(q)}, 
-V ys {q) = iU 2 [ Xzx {q) - Xxz(q)}- 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 3 




0.5 0.6 0.7 0.8 0.9 1.0 

U/U, 



Fig. 1. Interaction dependence of maximum eigenvalue for super- 
conductivity belonging to T-representation in Cav point group of 
the non-ccntrosymmetric system. 



gk a {ip(k)ip-(k) + d(k) • gk^+(k)]/2 are obtained with 
<p±(k)=tanh ^t/2E k+ ± tanh ^/2£ k _. 

In the centrosymmetric case, only V ss (q) and V xx (q) = 
V yy (q) remain, and the usual spin fluctuations in V ss (q) 
and V xx (q) contribute to the spin-singlet and spin-triplet 
pairing mechanisms, respectively. For non-zero Rashba 
field, the anomalous spin fluctuations contribute only 
to the spin-triplet pairing interactions V^^q) and the 
mixing interactions V sa (q) = — V as (q) for a,(3=x, y. 
We also note that the mixing interactions between spin- 
singlet and spin-triplet pairs are described by antisym- 
metric spin fluctuations, x yz {q) - Xzy(<l) and Xzx{<i) - 
Xxz(l), which relate with the Dzyaloshinski-Moriya in- 
teraction. 31 Therefore, the anomalous spin fluctuations 
are essential to form the spin-triplet pairs in the non- 
centrosymmetric system, if they are constructive for su- 
perconductivity. 

In the following, we show the results of numerical cal- 
culation. When the system is in the normal state, the gap 
equation reduces to an eigenvalue problem. When the 
maximum eigenvalue of a representation reaches unity, 
superconductivity belonging to this representation ap- 
pears. For actual calculation, we introduce a cutoff en- 
ergy u c =0.2ti for electrons forming the Cooper pair. Fur- 
thermore, we fix the superconducting transition tempera- 
ture at T c =0.Q2ti. In Fig. 1, it is shown that by increasing 
the on-site repulsion U toward an incommensurate mag- 
netic instability at C/ C =2.551ti, 18 every maximum eigen- 
value Ar of T-representation for the eigenvalue problem 
increases, and Xa 1 reaches unity first among all tetrag- 
onal group representations. This means that the super- 
conductivity of A\ representation in Ca v appears around 
the incommensurate magnetic instability. 

For the eigenvector of \a x =1 at the critical interaction 
constant Ua 1 , k-dependences of the corresponding nor- 
malized gap functions along Fermi-surfaces of e k _ and 
£k+ are shown in Figs. 2(a) and 2(b), respectively. Ig- 
noring all anomalous spin fluctuations with U=Ua 1 in 
the gap equation, the k-dependences of the gap functions 
change to Figs. 2(c) and 2(d), respectively, where Xa ± de- 
creases to 0.573 although the same interaction constant 
as Figs. 2(a) and 2(b) are used. Comparing these fig- 
ures, the singlet gap function ip(k) is not affected by the 
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Fig. 2. (a) and (b): k-dependence of gap functions, «/<(k) (solid), 
d x (]i) (dotted), and d y (k) (dashed) along Fermi-surfaccs of £k± 
for the full gap equation with the critical interaction constant 
satisfying A.a 1 =1. The insets show Fermi-surfaccs in the 1 st - 
quadrant of Brillouin zone, (c) and (d): k-dependence of gap 
functions along Fermi-surfaces of £k± for a gap equation with- 
out any anomalous spin fluctuations using the same interaction 
constant in (a) and (b). 




Fig. 3. k-dcpcndcncc of gap functions along Fermi-surface 

°f e k± a t several temperatures below T c . Inset shows the temper- 
ature dependences of spin-singlet and spin-triplet gap functions 
at momenta on Fermi-surfaccs giving the corresponding maxi- 
mum values of i/>(ki) (solid) and d(k2) ■ gk 2 (dashed), where kj 
are shown in insets of Fig. 2. 



anomalous spin fluctuations, as expected from eq. (0.10). 
Therefore, it is stabilized by the usual spin fluctuations 
with a momentum corresponding to the ordering wave 
vector at U c , which spans from one peak position of ip(\s) 
along a Fermi-line £k-=0 in Fig. 2(a) to other peak po- 
sition of |V>(k)| along a Fermi-line £k+=0 in Fig. 2(b). 18 
On the other hand, the magnitude of spin-triplet gap 
functions are enhanced by switching on the anomalous 
spin fluctuations, and becomes almost the same size as 
that of spin-singlet gap. Thus, the anomalous spin fluc- 
tuations are constructive for the superconductivity, and 
they are responsible for the stabilization of spin-triplet 
gap function. We note that the effect of anomalous spin 
fluctuations on superconductivity does not depend on ui c . 
Due to the nodal structures and similar amplitudes of 
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V>(k) and d(k), nodal structures for Ak± are also ex- 
pected. Considering that Ak± are the gap functions in 
the band picture, they reflect directly to thermal and dy- 
namical quantities in the superconducting state. For var- 
ious temperatures below T c =0.02ii, the k-dependences 
of Ak± along Fermi surfaces of £k± are shown in Fig. 3, 
where pairing interactions are fixed to those at T = T c . 
We stress that the k-dependences of gap functions re- 
flect fully the characteristic momentum dependences of 
spin fluctuations without assuming the common simple 
form of Ak± = A s ± A t |gk|. Recalling that the repre- 
sentation of the superconductivity is A\ in d v , both gap 
functions Ak± exhibit accidental line nodes. The posi- 
tions of nodes almost do not move with decreasing tem- 
perature, and the temperature dependence of V'(ki) and 
d(k2) • gk 2 is of common BCS type as shown in the inset. 
Therefore, the existence of the accidental line node in 
gap function explains T-dependences of NMR relaxation 
rate, thermal conductivity, and penetration depth in the 
superconducting state of CePt 3 Si and CeIrSi 3 . 8 ~ 10 ' 24 In 
order to check our scenario, it is desirable to observe the 
gap function in detail by the thermal conductivity exper- 
iment, in addition to the observation of anomalous spin 
fluctuations. 

Finally, we comment on the superconductivity in non- 
centrosymmetric correlated electron systems. As already 
mentioned, the spin-singlet and spin-triplet gap func- 
tions coexist in the non-centrosymmctric superconduc- 
tor. If only the usual antiferromagnetic spin fluctuations 
are enhanced, the spin-triplet gap function will be in- 
duced by the primary order parameter corresponding to 
the spin-singlet gap function, as shown in Figs. 2(c) and 
2(d). 16 However, in the case that the anomalous spin 
fluctuations are also enhanced, 18 the superconducting 
state is quite different from the former case, although 
it will not be so different from that of the usual uncon- 
ventional line-node superconductivity. With respect to 
its relation to magnetism, we note that unconventional 
superconductivity appears due to the spin fluctuations 
enhanced around the magnetic instability in a similar 
manner as in a centrosymmetric superconductor. On the 
other hand, the spin-triplet gap function develops due 
to the anomalous spin fluctuations characteristic of the 
non-centrosymmetric structure. 

In summary, we have studied the mechanism of super- 
conductivity in the non-centrosymmetric system. It has 
been shown that the spin-triplet gap function surviving 
the suppression of pairing due to the absence of cen- 
trosymmetry is developed by the anomalous spin fluctu- 
ations, which are enhanced around a magnetic instabil- 
ity. Reflecting the anisotropic momentum dependences 
of anomalous spin fluctuations, the gap functions have 
non-trivial Ai k-dependence with accidental line node. 
This pairing mechanism gives a reasonable explanation 
for the superconductivity in a non-centrosymmctric su- 
perconducting systems CePt 3 Si and CeTSi 3 , showing a 



huge upper critical field, a line-node gap structure, and 
an incommensurate magnetic structure. 
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